Entanglement between Collective Operators in the Linear Harmonic Chain 
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We investigate entanglement between collective operators of two blocks of oscillators in an infinite 
linear harmonic chain. These operators are defined as averages over local operators (individual 
oscillators) in the blocks. On the one hand, this approach of "physical blocks" meets realistic 
experimental conditions, where measurement apparatuses do not interact with single oscillators 
but rather with a whole bunch of them, i.e., where in contrast to usually studied "mathematical 
blocks" not every possible measurement is allowed. On the other, this formalism naturally allows 
the generalization to blocks which may consist of several non-contiguous regions. We quantify 
entanglement between the collective operators by a measure based on the Peres-Horodecki criterion 
and show how it can be extracted and transferred to two qubits. Entanglement between two blocks 
is found even in the case where none of the oscillators from one block is entangled with an oscillator 
from the other, showing genuine bipartite entanglement between collective operators. Allowing the 
blocks to consist of a periodic sequence of subblocks, we verify that entanglement scales at most with 
the total boundary region. We also apply the approach of collective operators to scalar quantum 
field theory. 

PACS numbers: 03.65.-w, 03.67.-a, 03.67.Mn 



I. INTRODUCTION 

Quantum entanglement is a physical phenomenon in 
which the quantum states of two or more systems can 
only be described with reference to each other, even 
though the individual systems may be spatially sepa- 
rated. This leads to correlations between observables 
of the systems that cannot be understood on the ba- 
sis of classical (local realistic) theories p). Its impor- 
tance today exceeds the realm of the foundations of quan- 
tum physics and entanglement has become an impor- 
tant physical resource, like energy, that allows perform- 
ing communication and computation tasks with efficiency 
which is not achievable classically 0]. 

In the near future we will certainly see more and more 
experiments on entanglement of increasing complexity. 
Moving to higher entangled systems or entangling more 
systems with each other, will eventually push the realm of 
quantum physics well into the macroscopic world. It will 
be therefore important to investigate under which con- 
ditions entanglement within or between "macroscopic" 
objects, each consisting of a sample containing a large 
number of the constituents, can arise. 

Recently, it was shown that macroscopic entanglement 
can arise "naturally" between constituents of various 
complex physical systems. Examples of such systems are 
chains of interacting spin systems 0-0- harmonic oscilla- 
tors 0,IH and quantum fields |{J. Entanglement can have 
an effect on the macroscopic properties of these systems 
000 and can be in principle extractable from them 
for quantum information processing 0, 0, ^] . 

With the aim of better understanding macroscopical 



entanglement we will investigate entanglement between 
collective operators in this paper. A simple and natural 
system is the ground state of a linear chain of harmonic 
oscillators furnished with harmonic nearest-neighbor in- 
teraction. The mathematical entanglement properties of 
this system were extensively investigated in [j, y, [ljj LL-J • 
Entanglement was computed in the form of logarithmic 
negativity for general bisections of the chain and for con- 
tiguous blocks of oscillators that do not comprise the 
whole chain. It was shown that the log-negativity typi- 
cally decreases exponentially with the separation of the 
groups and that the larger the groups, the larger the max- 
imal separation for which the log-negativity is non-zero 
0. It also was proven that an area law holds for har- 
monic lattice systems, stating that the amount of entan- 
glement between two complementary regions scales with 
their boundary |15| . 

In a real experimental situation, however, we are typ- 
ically not able to determine the complete mathemati- 
cal amount of entanglement (as measured, e.g., by log- 
negativity) which is non-zero even if two blocks share 
only one arbitrarily weak entangled pair of oscillators. 
Our measurement apparatuses normally cannot resolve 
single oscillators, but rather interact with a whole bunch 
of them in one way, potentially even in non- contiguous 
regions, thus measuring certain global properties. Here 
we will study entanglement between " physical blocks" of 
harmonic oscillators — existing only if there is entan- 
glement between the collective operators defined on the 
entire blocks function of their size, relative dis- 

tance and the coupling strength. Our aim is to quantify 
(experimentally accessible) entanglement between global 
properties of two groups of harmonic oscillators. Surpris- 
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ingly, we will see that such collective entanglement can 
be demonstrated even in the case where none of the oscil- 
lators from one block is entangled with an oscillator from 
the other block (i.e., it cannot be understood as a cumu- 
lative effect of entanglement between pairs of oscillators) , 
which is in agreement with p|. This shows the existence 
of bipartite entanglement between collective operators. 

Because of the area law |l5j the amount of entangle- 
ment is relatively small in the first instance. We suggest 
a way to overcome this problem by allowing the collective 
blocks to consist of a periodic sequence of subblocks. Then 
the total boundary region between them is increased and 
we verify that indeed a larger amount of entanglement is 
found for periodic blocks, where the entanglement scales 
at most with the total boundary region. We give an an- 
alytical approximation of this amount of entanglement 
and motivate how it can in principle be extracted from 

the chain H 03 OH [12 • 

Methodologically, we will quantify the entanglement 
between collective operators of two blocks of harmonic 
oscillators by using a measure for continuous variable sys- 
tems based on the Peres-Horodecki criterion 

EHG3 

fisj . The collective operators will be defined as sums over 
local operators for all single oscillators belonging to the 
block. The infinite harmonic chain is assumed to be in 
the ground state and since the blocks do not comprise 
the whole chain, they are in a mixed state. 



II. LINEAR HARMONIC CHAIN 

We investigate a linear harmonic chain, where each 
of the N oscillators is situated in a harmonic potential 
with frequency u> and each oscillator is coupled with its 
neighbors by a harmonic potential with the coupling fre- 
quency f2. The oscillators have mass m and their posi- 
tions and momenta are denoted as q i and p i , respectively. 
Assuming periodic boun dary conditions (<?jv+i = 5i)> the 
Hamiltonian thus reads j20j 



H = 



N 

E 



2 m 



m lo Qj mQ, 2 {q 3 - q^) 



(1) 



We canonically go to dimensionless variables: qj = C qj 

and p 3 = pj/C, where C = y/mto(l + 2 Q 2 /lo 2 ) 1/2 0- 
By this means the Hamiltonian becomes 



p N 



(2) 



with the abbreviations a = 2$} 2 /(2fl 2 + uj 2 ) and E = 
y/2 fl 2 + lo 2 . The coupling constant is restricted to val- 
ues < a < 1, where a — > in the weak coupling 
limit (Q/lo — ► 0) and a — > 1 in the strong coupling limit 
{Vl/lo — > oo). 

In the language of second quantization the positions 
and momenta are converted into operators {qj — > qj, 



Pj — > pj) and are expanded into modes of their anni- 
hilation and creation operators, a and a\ respectively: 



N-l 



^ = -^E^7^[«^)e ie ^" + H.c.], (3) 



N-l 



[a(9 k )e ie »J -H.c.]. (4) 



Here 9 k = 2-rrk/N (with k = 0, 1, ...,N-1) is the dimcn- 



sionless pseudo- momentum and v(9k) = \/l — a cos 9 k 
is the dispersion relation. The annihilation and cre- 
ation operators fulfil the well known commutation re- 
lation [a(9k),al(9k')] = 5 kk >, since [qupj] = i<% has to 
be guaranteed. The ground state (vacuum), denoted as 
|0), is defined by a(9 k ) |0) = holding for all 9 k . The 
two-point vacuum correlation functions 



9\ 



I <l i 1j 



= (0\PiPj |0) = (PiPj) , 



(5) 
(6) 



are given by gi = (2 N) 1 J2 k =o 1 ' 1 (^k)cos(l9 k ) and 
hi - (2N)- 1 J2^ 1 ^9 k )cos(l9 k ), where I = \i - j\. 
In the limit of an infinite chain (TV — > oo) — which 
we will study below — and for I < N/2 they can 
be expressed in terms of the hypergeometric function 
2 F 1 0: 9l = (z l /2n)( l -y 2 ) 2*1(1/2,1+ 1/2,/ + l,z 2 ), 
hi = (/x z l /2)( 1 - 3 / 2 ) 2*x (-1/2, / - 1/2, 1 + 1, z 2 ), where 
z = (1 - Vl - a 2 ) I a and /j, = + z 2 . 



III. DEFINING COLLECTIVE OPERATORS 

In the following, we are interested in entanglement 
between two "physical blocks" of oscillators, where the 
blocks are represented by a specific form of collective op- 
erators which are normalized sums of individual oper- 
ators. By means of such a formalism we seek to fulfil 
experimental conditions and constraints, since finite ex- 
perimental resolution implies naturally the measurement 
of, e.g., the average momentum of a bunch of oscillators 
rather than the momentum of only one. On the other 
hand, this formalism can easily take account of blocks 
that consist of non- contiguous regions, leading to inter- 
esting results which will be shown below. We want to 
point out that this convention of the term block is not the 
same as it is normally used in the previous literature. In 
contrast to the latter, for which one allows any possible 
measurement, our simulation of realizable experiments 
already lacks some information due to the averaging. 

Let us now consider two non-overlapping blocks of os- 
cillators, A and B, within the closed harmonic chain in 
its ground state, where each block contains n oscillators. 
The blocks are separated by d > oscillators (Fig. ^| . 
We assume n, d <C N and N — > oo for the numerical 
calculations of the two-point correlation functions. 
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FIG. 1: Two blocks of a harmonic chain A and B. Each 
block consists of n oscillators and the blocks are separated by 
d oscillators. 

By a Fourier transform we map the n oscillators of each 
block onto n ("orthogonal") frequency- dependent collec- 
tive operators 

with the frequencies A; — 0, n — 1, and analogously for 
block B. The commutator of the collective position and 
momentum operators is 

[Qf,P^] = r8 w . (9) 

This means that collective operators for different frequen- 
cies k 7^ k' commute. For different blocks the commuta- 
tor vanishes: [Qa\Pb' ] ] = 0. 

If the individual positions and momenta of all oscilla- 
tors are written into a vector 

x = (q 1 ,p 1 ,q 2 ,p2, ...,q N ,p N ) T , (10) 

then there holds the commutation relation 

[xi,Xj] = idij (11) 

with $1 the n-fold direct sum of 2x2 symplectic matrices: 

;)■ (i2) 

3=1 V 7 

A matrix S transforms x into a vector of collective (and 
uninvolvcd individual) oscillators: 

X ee Sx = ({Qf, Pf } k , {Q« Ff>} k , {qj^hf. 

(13) 

Here {Qf,Pf} k = (Qf, Pf , Q^^t de- 
notes all collective oscillators of block A and analogously 
for block B, whereas {qj,Pj}j denotes the 2 (N — 2n) 
position and momentum entries of those N — 2 n oscil- 
lators which are not part of one of the two blocks. The 
matrix S corresponds to a Gaussian operation J^. It 
has determinant det S = 1 and preserves the symplectic 
structure 

fJ = S T r&S, (14) 



and hence 

[jt i ,jt j ]=isi i 3 (is) 

for all in particular verifying This means that the 
Gaussianness of the ground state of the harmonic chain 
(i.e., the fact that the state is completely characterized by 
its first and second moments, see below) was preserved by 
the (Fourier) transformation to the frequency-dependent 
collective operators. 

IV. QUANTIFYING ENTANGLEMENT 
BETWEEN COLLECTIVE OPERATORS 

In reality, we are typically not capable of single parti- 
cle resolution measurements and only of measuring the 
collective operators with one frequency, namely k = 0. 
Note that in general the correlations of higher-frequency 
collective operators, e.g., ((Q^) 2 ) or (QaQb) with 
k =/= 0, are not real numbers. Therefore, as a natural 
choice, we denote as the collective operators 

Q A ^Qf = 4=£&, ( 16 ) 

^Pf = ^£^, (17) 

and analogously for block B. It seems to be a very nat- 
ural situation that the experimenter only has access to 
these collective properties and we are interested in the 
amount of (physical) entanglement one can extract from 
the system if only the collective observables Qa,b and 
Pa,b are measured. 

Reference [Tsj derives a separability criterion which is 
based on the Peres-Horodecki criterion [H E3 and the 
fact that — in the continuous variables case — the partial 
transposition allows a geometric interpretation as mirror 
reflection in phase space. Following largely the notation 
in the original paper, we introduce the vector 

i=(Q A ,P A ,Q B ,P B ) (18) 

of collective operators. The commutation relations 
have the compact form [£ a , = i K a p with K = 

2 =1 1 ^ jj. The separability criterion bases on the co- 
variance matrix (of first and second moments) 

V a0 = ^( + Aip A£ a ) , (19) 

where A£ Q = £ a — ( £ a ) with (£ a ) = in our case (state 
around the origin of phase space) . 

The covariance matrix V is real (which would not 
be the case for higher-frequency collective operators) 
and symmetric: (QaQb) = {QbQa) and (PaPb) = 
(PbPa), coming from the fact that the two-point cor- 
relation functions and © only depend on the 
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absolute value of the position index difference. On 
the other hand, using © and we verify that 

(QiPj) = i(2-/V)~ 1 Efclro lex P[ i ^(* ~ 3)] and (PjQi) = 
— i (2N)^ 1 J2k=o ex P[i^fe(j ~ *)]■ For i ^= j both summa- 
tions vanish (6k = 2nk/N and i,j integer) and for i = j 
they are the same but with opposite sign. Thus, in all 
cases (qipj) = — (pj qi ). These symmetries also hold for 
the collective operators and hence we obtain 



d = 



(a) 



(b) 



G Gab 

OHO Hab 

Gab G 

H AB H 



(20) 



The matrix elements are 



= = <<&> = - ££$0-4, (21) 

j£A ieA 

H={P%) = {P%) = \Y,T, h \o-i\^ ( 22 ) 

Gab = (QaQb) = ~^2^2g\j-i\ , (23) 



3'eA ieB 

\ \ ^ 



(24) 



jEA iG-B 



To quantify entanglement between two collective blocks 
we use the degree of entanglement e, given by the ab- 
solute sum of the negative eigenvalues of the partially 
transposed density operator: e = Tr|p Ts | — 1, i.e., by 
measuring how much the mirror reflected state fails to 
be positive definite. This measu re (neg ativity) is based 
on the Peres- Horodecki criterion Eo, QjJ an d was shown 
to be an entanglement monotone l22l 231 . For covariance 
matrices of the form 120(1 it reads [ljj 



max 



(SiS 2 )o 
SiS 2 



1 



(25) 



where S± = G — \Gab\ and S2 = H — \Hab\- In general, 
the numerator is defined by the square of the Heisenberg 
uncertainty relation 



(<VMn = [\\{[Q A ,B,Pam]\ 



(26) 



with (<5i <52)o = 1/4 due to (J5J). We note that £ is a degree 
of entanglement (in the sense of necessity and sufficiency) 
only for Gaussian states which are completely character- 
ized by their first and second moments, as for example 
the ground state of the harmonic chain we are studying. 
However, we left out the higher-frequency collective op- 
erators (and all the oscillators which are not part of the 
blocks) and therefore, the entanglement e has to be un- 
derstood as the Gaussian part of the amount of entangle- 
ment which exists between (and can be extracted from) 
the two blocks when only the collective properties Qa,b 
and Pa,b, as defined in l|ltj|) and (|17|l . are accessible. 




2 4 6 8 10 




FIG. 2: Degree of collective entanglement e for two blocks 
of oscillators as a function of their size n. (a) The blocks 
are neighboring (d = 0) and entanglement exists for all n 
and coupling strengths a. Plotted are a — 0.99 (diamonds), 
a = 0.9 (squares) and a = 0.5 (triangles), (b) The same for 
two blocks which are separated by one oscillator (d = 1). The 
two blocks are unentangled for n = 1 but can be entangled, 
if one increases the block size (n > 1), although non of the 
individual pairs between the blocks is entangled. 



There also exists an entanglement witness in form of 
a separability criterion based on variances, where A = 
((Qa-Qb) 2 ) + ((Pa+Pb) 2 ) =2(G-Gab+H+H ab ) < 
2 is a sufficient condition for the state to be entangled 
[24|. We note that the above negativity measure H25|) 
is "stronger" than this witness in the whole parameter 
space (a,n). In particular, there are cases where s > 
although A > 2. This is in agreement with the finding 
that the variance criterion is weaker than a generalized 
negativity criterion |25| . 

We further note that the amount of entanglement 
(|25|l is invariant under a change of potential redefini- 
tions of the collective operators, e.g., Qa = J2jeA^i 
or Qa = (l/fi) Y^jeA as then the modified scaling 
in the correlations (G, Gab, H, and Hab) is exactly 
compensated by the modified scaling of the Heisenberg 
uncertainty in the numerator. 

Figure |21 shows the results for d = and d = 1. In the 
first case — if the blocks are neighboring — there exists 
entanglement for all possible coupling strengths a and 
block sizes n. In the latter case — if there is one oscil- 
lator between the blocks — due to the strongly decaying 
correlation functions g an h there is no entanglement be- 
tween two single oscillators (n = 1), but entanglement for 
larger blocks (up to n = 4, depending on a). The state- 
ment that entanglement can emerge by going to larger 
blocks was also found in Q|. But there the blocks were 
abstract objects, containing all the information of their 
constituents. In the case of collective operators, however, 
increasing the block size (averaging over more oscillators) 
is also connected with a loss of information. In spite of 
this loss and the mixedness of the state, two blocks can 
be entangled, although non of the individual pairs be- 
tween the blocks is entangled — indicating true bipartite 
entanglement between collective operators. For d > 2, 
however, no entanglement can be found anymore. 

These results are in agreement with the general state- 
ment that entanglement between a region and its com- 
plement scales with the size of the boundary |l5j |. In 
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FIG. 3: Two periodic blocks of a harmonic chain A and B. 
Each block can consist of m subblocks with s oscillators each, 
separated by d oscillators. In the picture d= 1, m — 2, s = 3 
and the number of oscillators per block is n = ms = 6. 



the present case of two blocks in a one-dimensional chain 
(Fig. ^| the boundary is constant and as the blocks are 
made larger, the entanglement decreases since it is dis- 
tributed over more and more oscillators. We therefore 
propose to increase the number of boundaries by consid- 
ering two non-overlapping blocks, where we allow a peri- 
odic continuation of the situation above, i.e. a sequence 
of m > 1 subblocks, separated by d oscillators and each 
consisting of s > 1 oscillators, where ms = n (Fig. [3J. 

The degree of entanglement between two periodic 
blocks of non-separated (d = 0) one-particle subblocks 
(s = 1) is larger for stronger coupling constant a and 
grows with the overall number of oscillators n (Fig. 0^,). 
For given a and n and no separation between the sub- 
blocks (d — 0) the entanglement is larger for the case of 
small subblocks, as then there are many of them, causing 
a large total boundary (Fig. @J>). Entanglement can be 
even found for larger separation (d = 1,2) with a more 
complicated dependence on the size s of the subblocks. 
There is a trade-off between having a large number of 
boundaries and the fact that one should have large sub- 
blocks as individual separated oscillators are not entan- 
gled (Fig.^J;,d). For d > 3 no entanglement can be found 
anymore. (In a realistic experimental situation, where 
the separation d is not sharply defined, e.g., where there 
are weighted contributions for d — 0, 1, d max , entan- 
glement can persist even for d max > 3, depending on the 
weighting factors.) 

For the sake of completeness we give a rough approxi- 
mation of the entanglement between two periodic blocks. 
Let us assume that the subblocks are directly neighbored, 
d = 0. Furthermore, we consider couplings a such that 
we may neglect higher than next neighbor correlations 
(a < 0.5), i.e., we only take into account go, gi, ho and 
hi. The correlations read 



G 



2m(g-l) 
~ .9o H 9i , (27) 



Gab = ~ 9\i-i\ ~^( 2 m - !) 9i 



(28) 



j£A i£B 



and analogously for H and Hab- The first equation re- 
flects that there are n self-correlations and m (s— 1) near- 
est neighbor pairs (which are counted twice) within one 
block, i.e., s — 1 pairs per subblock. The second equa- 
tion represents the fact that there are 2 m — 1 boundaries 
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FIG. 4: Degree of collective entanglement e for two periodic 
blocks of oscillators as a function of their total size n. (a) 
Neighboring one-particle subblocks (d = 0, s = 1). Entan- 
glement monotonically increases with n and becomes larger 
as the coupling strength a increases. Plotted are a — 0.99 
(diamonds), a = 0.9 (squares) and a — 0.5 (triangles), (b) 
The coupling is fixed to a = 0.99 for this and the subse- 
quent figures. There is no separation, d = 0. Plotted are 
the cases s = 1, 2, 5. For fixed n the entanglement is more or 
less proportional to the number of boundaries, i.e., inversely 
proportional to the subblock size s. (c) and (d) correspond 
to the cases d — 1 and d — 2, respectively. The dependence 
on the size of the subblocks is more complicated as there is 
a trade-off between having a large number of boundaries (i.e. 
small s) and the fact that one should have large subblocks as 
individual separated oscillators are not entangled. 



where blocks A and B meet. Using s — n/m, the entan- 
glement (|25f) becomes (note that g\ > and hi < 0) 



1 



4\ go + (2-±*=±)g 1 ] [/i + (2-i)/n 



(29) 



For given n this approximation obviously increases with 
the total number of boundaries, m. It can be consid- 
ered as an estimate for a situation like in Fig. 0Jd, if a 
smaller coupling is used such that the neglect of higher 
correlations becomes justified. 

We close this section by annotating that the entan- 
glement l|25ll between collective blocks of oscillators - 
being the Gaussian part — can in principle (for suffi- 
cient control of the block separation d) be transferred 
to two remote qubits via a Jaynes-Cummings type in- 
teraction 0, IToL IT2l| . For the interaction with periodic 
blocks " gratings" have to be employed in the experimen- 
tal setup. The interaction Hamiltonian is of the form 



(e 



"l 



e 



(30) 



where u>i is the Rabi frequency and = (o - ")^ = |e)j i(g\ 
is the bosonic operator (with \g)i and |e)j the ground and 
the excited state) of the i-th qubit (i — 1,2). 
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V. COLLECTIVE OPERATORS FOR SCALAR 
QUANTUM FIELDS 

The continuum limit of the linear harmonic chain is the 
(l+l)-dimensional Klein-Gordon field 4>(x,t) with the 
canonical momentum field Tr(x,t) = <p(x,t). It satisfies 
the Klein-Gordon equation (in natural units K = c = 1) 
with mass m 



tered at position xq. The commutator is 



+ 171 = 0. 



(31) 



With the canonical quantization procedure (f> and tt be- 
come operators satisfying the non-trivial commutation 
relation [<f)(x, t), -k(x', t)] = iS(x — x'). The field operator 
can be expanded into a Fourier integral over elementary 
plane wave solutions |26| 



4>{x,t) 



dk 



V4 7T LO k 



[a(k)e ikx -' luJkt + H.c], (32) 



tt(x, t) = — i 



dk LUk 



a(k) e 



i k x— i lj/c t 



H.c], (33) 



where fc is the wave number and u>k = +y/k 2 + m 2 is the 
dispersion relation. The annihilation and creation oper- 
ators fulfil [a(k),al(k')] = 5(k - k'). We write the field 
operator as a sum of two contributions <\> — + 4>^~~\ 
where (f>^ (<p is the contribution with positive (nega- 
tive) frequency. Thus, <p + ) corresponds to the term with 
the annihilation operator in 1)32(1. The vacuum correla- 
tion function is given by the (equal-time) commutator of 
the positive and the negative frequency part: 

(0| j>(x, t) 4>(y, t) |0) = [$M (x, i), (y, t)] . (34) 

It is a peculiarity of the idealization of quantum field 
theory that for x = y this propagator diverges in the 
ground state: 



(j) 2 {x,t) |0) — ► oo . 



(35) 



The same is true for (0| tt 2 (x, t) |0) and hence we cannot 
easily build an entanglement measure like for the har- 
monic chain, since the analogs of the two-point correla- 
tion functions go and ho, © and 10, are divergent now. 
Automatically, we are motivated to study the more phys- 
ical situation and consider extended space-time regions, 
which means that we should integrate the field (and con- 
jugate momentum) over some spatial area. We define the 
collective operators 



1 



L/2 
L J-L/2 

1 r L l 2 



(x + Xo , t) dx , (36) 



R L (xo,t) = —7= / tt{x + x ,t)dx , (37) 

VL J-L/2 

Therefore, &L{xo,t) and U.L(xo,t) are equal-time oper- 
ators which are spatially averaged over a length L, cen- 



L/2 r L/2 



L/2 J-L/2 



id^\x-y)=i, 
(38) 

which is in complete analogy to If $l and IT^ 

correspond to separated regions without overlap, i.e., 
\x - 2/0 1 > L, then of course [$ L (xo,t),fl L (yo,t)] = 0. 
The spatial integration in l(36|l and 1(37(1 can be carried 
out analytically: 



®L(x ,t) 



fl L (x ,t) 




dk fkL 
- — — — sin —— 

k JuJk \ 2 



i k XQ — iujk t 



H.c], 



<t:L 

x [a(k)e ikx °- iukt 



dk JuJk fkL 

sm 

k V 2 



H.c.]. 



(39) 



(40) 



The final step is to calculate the propagators of the field 
and the conjugate momentum. We find 



$L(x ,t) $ L (yo,t)\0) 



(41) 



I r dk . 2 (kL\ , 

/ ■ sin coslfc r) . 

nLj_ QO k 2 Vk 2 + m 2 \ 2 1 



D AL (r) = (0\IL L (xo,t)U L {y ,t)\0) 



(42) 



1 

7rL 



dk \/k 2 + m 2 
k 2 



kL 



sin — - I cos(fcr) , 



with r = \xq — yo I the distance between the centers of 
the two regions, reflecting the spatial symmetry. Thus 
Dj, L (0) and L (0) are the analogs of (Q 2 A B ) and 

(P\ B ) (iniro-block correlations within the same block), 
respectively, whereas L (r > L) and L (r > L) cor- 
respond to (QaQb) and (PaPb) (inter-block correla- 
tions between separated blocks). 

The expressions ((4111 and 1(42(1 are finite, especially for 
r = 0. Mathematically, the integration over a finite spa- 
tial region L corresponds to a cutoff, which removes the 
divergence we faced in (135(1 . However, the expressions 
are ill defined for L — > 0. 

Applying the entanglement measure 1(25(1 with G = 
D $L (0), H = D^O), G AB = D $L (r), and H AB = 
Djj L (r) does not indicate entanglement for any choice of 
L and r > L. The same is true for the generalized case of 
blocks consisting of periodic subregions of space, showing 
an inherent difference between the harmonic chain and its 
continuum limit. We believe this is due to the fact, that 
any spatial integration immediately corresponds to an 
infinitely large block in the discrete harmonic chain and 
that the information loss (compared to the mathematical 
indeed existing exponentially small entanglement |6j) due 
to the collective operators already is too large. Nonethe- 
less, defining collective operators like in (|S5|l and (l3*7|) 
and use of the measure ((25(1 may reveal entanglement 
between spatially separated regions for other quantum 
field states, which is the subject of future research. 
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VI. CONCLUSION 

Our results have importance for investigating the con- 
ditions under which entanglement can be detected by 
measuring collective observables of blocks consisting of a 
large number of harmonic oscillators. This has relevance 
for schemes of extracting entanglement where the probe 
particles normally interact with whole (periodic) groups 
of oscillators rather than single oscillators. The results 
are also relevant for the transition from the quantum to 
the classical domain as they suggest that entanglement 
between collective operators (global properties) may per- 
sist even in the limit of a large number of particles. It is 
obvious that our approach of collective observables can be 
extended to more dimensions. Furthermore, we demon- 



strated its potential application to scalar quantum field 
theory. 
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